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ABSTRACT 

Large scale surveys of the prominent members of the Local Group have provided compelling evidence 
for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought 
to be the debris from ancient and on-going accretion events. In this paper, we compare two extant 
surveys of the M31-M33 subgroup of galaxies; the Pan- Andromeda Archaeological Survey (PAndAS) 
of the stellar structure, and a combination of observations of the H I gaseous content, detected at 
21cm. Our key finding is a marked lack of spatial correlation between these two components on all 
scales, with only a few potential overlaps between stars and gas. The paucity of spatial correlation 
significantly restricts the analysis of kinematic correlations, although there does appear to the H I 
kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results 
demonstrate that that different processes must significantly influence the dynamical evolution of the 
stellar and H I components of substructures, such as ram pressure driving gas away from a purely 
gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will 
provide a determination of the properties of the gaseous halo of M31 and M33. 

Subject headings: galaxies: abundances, galaxies: dwarf, galaxies: kinematics and dynamics, Local 
Group, dark matter 



1. INTRODUCTION 

In ACDM cosmologies, galactic halos are built up over 
time through the accretion and cannibalisation of smaller 
mass systems. Given the complexity of the dynami- 
cal interactions, these are studied in detail using com- 
puter simulations (e.g. Ibata fe Lewis|l998 Fardal et al. 
2008 1 , although, due to the difficulties m modelling bary 



omc physics, these typically only consider the evolution 
of dark matter. However, these predict that the outer 
region of galactic halos should be dominated by exten- 
sive tidal streams, whereas the shorter dynamical time 
scales of the inner halo result in com plete destruction 
and a smooth stellar distribution (e. g. [B ullock fe John 



ston|2005| |Cooper et al]|2010[ |Font et al.||201l| ) 



Trie progenitors of substructure could possess a mix of 
stars and gas, generally more concentrated than the dark 
matter halos in which they reside. The precise composi- 
tion of any in-falling halo will depend upon its mass and 
evolutionary history, and we would expect dwarf galax- 
ies to be stripped of their gas due to interactions within 
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the Local Group, whereas larger galaxies will hold-on to 
its gaseous component. While stars are effectively col- 
lisionlcss, more complex internal physics influences the 
evolution of gas, including internal shocking, cooling and 
collapse, star formation and its associated feedback, as 
well as ram-pressure strippi ng due to the presence of hot 



galactic halo gas (e.g. see Moore fc Davis 1994 Mu- 



raka mi fc Babul] [19991 |Bekki||2008| |Mastrop etro et al 
|2009[ ). This is quite apparent in the two major accre- 
tion events within the Milky Way halo, with the body 
and stream of the Sagittarius Dwarf Galaxy possessing 
no gas, having potenti ally lost it during its initial interac- 
tion with the Gala xy (Burton fc Lockman||1999 but see 
also Putman et al. 2004[ ), while the Magellanic Stream 
appears to be completely devoid of stars, consisting of 
stripped gas, although the question of whether this is the 
result of tidal or ram pressure forces has yet to be decided 
(IStanimirovic et al. ,2008; Be sla et al.|2010[ | Diaz fc Bekki 
201 1\ . However, while we should expect dmerent distri- 
butions for various disrupted components in M31 and 
M33, the complexities of ram-pressure stripping, which 
may req uire full magneto-hydro dynamic approaches to 
simulate ( Ruszkowski et al.|2012 ) , makes drawing robust 
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Fig. 1. — A schematic illustration of the prominent substructure apparent in Red Giant Branch (RGB) stars; this encompasses a broad 
cut in metallicity and nonlinearly scaled to enhance the stellar substructure. The thick solid line represents the PAndAS footprint, with 
the thinner solid lines resenting the main stellar substructure; as well as structure in the halo of M31, there is also tidal debris associated 
with M33 and NGC147. A number of substructures have been presented in previous publications (see Section |2.1.2[ l, but several additional 
features are apparent in t his full map, including t he Eastern Cloud and the broad swath of stars associated with Stream D (s ee also the 
PAndAS map presented in |Richardson et al.|2011[ | . The dashed line represents an over-density of globular clusters identified in [Mackey et| 
|al.| ( |2010[) . Closer in to the disk ot 1V131, tne dotted lines represents the inner shells and loops thought to be wraps of the Giant Stellar 
Stream. The large dot-dashed circle is at a radius of 150 kpc from the centre of M31, and the smaller being 50 kpc from the centre of M33. 



conclusions difficult. 

We are now in an era where large scale surveys are pro- 
viding a global picture of the formati on and evolution 
of galaxies over cosmic time (e.g. see Abraham fc van] 
den Bergh 2001). Given their distances, unraveling the 
fine details of galactic evolution is below the resolution 
achievable with the vast majority of galaxies in the Uni- 
verse. This is unfortunate, as it is on small scales, with 
the structure within the halos of large galaxies, that have 
presented the major challenge to the prevaili ng ACDM 



al. 1999 



jucki 



Moore et al."1|1999| ) 



paradigm (e.g. the missing satellite problem; Klypin et 

T 



y, the large members of the Local Group, the 
Milky Way, the Andromeda (M31) and Triangulum 
(M33) galaxies, are close enough for such hne scale detail 
to be resolvable, although their immense angular scale 
presents a significant challenge to the observability of 
their extensive halos. To this end, the last decade has 
seen the advent of large, detailed surveys of both the stel- 
lar and gaseous components of Local Group members, for 
the first time providing a panoramic view of galaxy for- 
mation in action; for the Milky Way, this includes sur- 



2011), SEGUE 



veys such a s RAVE flSiebert et al. 
Jong et al.||2010[ ), GAtiti ( |McGlure-Grimt'hs et alpot 
and others. 

In this paper, we present a detailed comparison of the 



stellar and gaseous matter in the halos of our nearest 
large companions in the Local Group, M31 and M33, us- 
ing a new map of stellar structure obtained through the 
Pan-Andromeda Archaeological Survey (PAndAS) and 
the most detailed maps of H I, obtained during several 
telescope campaigns. In Section [2j we outline the obser- 
vational campaigns that resultedin the data presented, 
with a summary of their key scientific results. In Sec- 
tion [3] we compare the spatial and kinematic distribu- 
tions of gaseous and stellar substructure. We present 
the interpretation of these correlations in Section [4] and 
our conclusions in Section [5l 

2. OBSERVATIONS & PROPERTIES 

Some galaxies of the Local Group hav e been known 
since pre-history (e.g. see Section 3.1 of Geehan et al. 
2006 and have played a pivotal role in our explorations 
ot the workings of the Universe. Reproducing a detailed 
account of this history is beyond the scope of this arti- 
cle, and the interested reader is directed to recent reviews 

3 Given its proximity, M31 is one of the few galaxies visible 
to the naked eyes. Its existence was first documented by the Per- 
sian astronomer Abd al-Rahman Al-Sufi in his treatise on stellar 
astronomy titled Kitab al-Kawatib al-Thabit al-Musawwar (Book 
on the Constellations of the Fixed Stars), published in AD 964, 
where he both identified its position in the sky and summarized 
his observations. 
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Fig. 2. — The distribution of H I gas (red contours), drawn from the survey of |Braun fc Thilker| ( fe004|l, overlaid on a surface density map 
of stars identified as being on the RGB within M31 and M33 in the entire PAndAS footprint (sec McUonnachic ct al. 2009 Richardson ct 
|al.|201l| . The logarithmic contours are drawn at integrated column densities of H I, with N^j/ = 10 iY to lO^" cm z in steps of 0.5 dex, 
and the annotations indicate the heliocentric radial velocity of H I features. The stellar map again represents a broad range of metallicitics 
and has a non-linear stretch to accentuate the substructure. 



e.g. |Tolstoy et aL 



Connachic ct al. 2 U12ap 



Wyse||2010| |Tolstoy||2011| |Mc- 
In the following we will focus 
upon the recent optical and radio surveys of M31 and 
M33, with a review of the observational programs and 
key scientific discoveries to date. 

2.1. Stellar Observations 
2.1.1. Observational Program 

The Pan- Andromeda Archaeological Survey (PAndAS) 
is a survey of 400 square degrees, covering the halo of 
M31 out to a distance of 150 kpc, and M33 to a cor- 
responding distance of 50 kpc, undertaken as a large 
program on MegaCarrj^J mounted on the 3.6-m Canada- 
France-Hawaii Telescope (CFHT). Integrations were suf- 
ficient to achieve photometric limits g = 25.5 and i — 
24.5 at a S/N ~ 10, reaching several magnitudes below 
the tip of the Red Giant Branc h (RGB) at the dist ance of 

The 
2011, 



2011 



M31/M33 (D- 780 - 900 kpc; |Conn et al 
photometric data-taking was completed in early 
with the first published map, covering roughly half of 
the to tal observed area, presented in McConnachi e et aT] 
fl2009|), with the (almost ) entire dataset first presented in 
Richardson et al. (2011). The final map and high level 



4 cf ht .hawaii . edu/Instruments/Imaging/Megacam 

5 In the following, we adopt the distances £>M3i = 779 kpc and 
D mh = 820 kpc, eac h with an uncertainty of ± 20 kpc, as derived 
by jConn et al.| l |2012| . 



data products will b e made publicly available in a for th- 
commg publication ( |McConnachie et al.|[2012a| |2013[) 

This study focuses upon Red Giant Branch (RGB) 
stars at the distance of M31 and M33, selected with cuts 



in color and magnitude (see Ibata et aL||2001 2007 Mc 



Connachie et al. 2009). Figure [T] presents the distribu- 



tion ot RGB stars, overlain by a schematic map of the 
prominent stellar substructure; note that a non-linear 
scaling has been applied to the RGB density, to bring out 
faint substructure, and contains a broad swath of metal- 
licities to reveal metal-poor (—3.0 < [Fe/H] < —1.7), 
intermediate (—1-7 < [Fe/H] < —0.7) and metal-rich 



0.7 < [Fe/H] < 0.0) substructure (cf. 



Ibata et al. 



2007] . The thick solid line is the entire PAndAS foot- 
prmt, with stellar substructure as labelled thin solid 
lines. The dashed curve represents a significant over- 
de nsity of globular clu sters identified as the NW Group 



by Mackey et al. (2010). The large dot-dashed line corre- 
sponds to a circle ot radius of 150 kpc from the center of 
M31, whereas the smaller dot-dashed circle represents a 
distance of 50 kpc from the center of M33. It is very ap- 
parent that accompanying M31 and M33 is a wealth of 
substructure consisting of extensive streams and dwarf 
gala xies [t hese wil l be discussed in more detail in Sec 
tion 



2.1.2 and see Richardson et al. 



The stellar catalog derived from 



pOTT )]. 
T^ndAS 



built 



upon earlier observat ions with the CFHT/MegaCam fc 



CHFT/CFH12k (see McConnachie et al. 2003 Ibata et 
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al. 2007). In parallel, a number of fields, targeting 
prominent substructure and dwarfs, were targeted with 
DEIMOS (Fa ber et al.1[2003| on the 10m Keck-II Tele - 
scope (e.g. |Chapman et al. 2008 Collins et al. 2011). 
With a moderate resolution [K ~ 6000), observations of 
60-90 mins around the prominent CaT absorption lines 
~ 8600A resulted in a S/N~ 5 for targets of i ~ 21, with 
a corresponding velocity resolution of ^ 5 — 10 km s _1 ; 
these will be described in more detail in Section I3T2I 

2.1.2. Stellar Properties 

Detailed analysis of the stellar subst ruct ure in PAndAS 
will be prese nted in Ibata et al. (2012) and McConnachie 
et al. ( |2013 ), so here we focus on the key features. 

There is a wealth of substructure apparent in the outer 
disk/inner halo regions of M31, revealed in the earlier 
Wide Field Camera o n the 2.5-m Isaac Newton Telescope 



as well as a new dwarf galaxy, Cass II (also known as 
And XXX), which is a potential satellite of the 147/185 
subgroup, and will b e discussed in mor e detail in a forth- 
coming publication ( Irwin et al.||2012 ). 



2.2. HI Observations 
2.2.1. Observational Program 

The ex tended environment of M 31 and M33 was ob- 
served by |Braun fc Thilker| p00T[) using the Westerbork 
Synthesis Radio Telescope (WSRT) as fourteen single 
dish telescopes. A region of 60 x 30° in RAxDec was 
imaged in the H I emission line at an effective resolu- 
tion of 49 arcmin with an RMS sensitivity corresponding 
to a column density of 5 x 10 16 cm~ 2 over ~ 30 km s _1 . 
This unprecedented sensitivity permitted detection of 
H I in emission from column densities that have previ- 
ously only been probed by Lya absorption toward back- 
ground QSOs. This low resolution survey has been sup- 
plemented by higher resolution total power observations 
made with the Green Bank Telescope (GBT). A region 
covering 7° x 7° centered on M31 and 5° x 5° centered on 
M33 was observed with GBT during 2002 October with 
multiple perpendicularly scanned coverages yielding im- 
ages with angular resolution as high as 9 arcmin. The 
M3 1 data have been presented previously by |Thilker et| 
Accompanying the GSS is a series of streams, labelled [iTl (12004} , while the M33 data were obtained wit 



surve y of the galaxy J phata et al.||2001 Richardson et al. 

[Tl), with the most significant struc 



2008[ & see Figure 

ture being an extensive stellar stream, the Giant Stellar 
Stream (GSS), in the halo of M31. This is seen as a 
metal-rich substructure wrapped within a more metal- 



M31 in the north-east (e.g. 


Ferguson et al. 2002 


Ibata et 


al. 


2007 


Fardal et al. 


2012 


1. 



B, C and D, which lie perpendicular to the Giant Stream 
(Ibata et al. 2007), plus Stream A, a stellar over-density 
~ 125 kpc (m projection) from M31. Further features 
include the Sou th West (SW) Cloud and the North West 
(NW ) Stream ( |McConnachie et al.|[2009[ |Mackey et~aL 
2010|) . As revea led in Figure I] and as shown by Richard- 
son et al. (2011 ), this latter feature appears to loop back 
towards M31 as the North East (NE) Stream, and inter- 
cepts the A nd XXVII dwarf galaxy; following |Carlberg] 
et al. (2011 ), we identify this entire feature as the North 
West Stream, and label the two components of it as NW1 
and NW2. Finally, also apparent in this image is another 
significant over-density of stars, named the Eastern (E) 
Cloud, located a distance of ~ 125 kpc from M31. 

As well as the extensive substructure, the PAndAS ob- 
servations revealed the p resence of a smoot h stellar halo 

see also the 



th the 

same setup and have been reduced in a similar fashion. 
These moderate resolution data are supplemented with 
interferometric mosaic observati ons consisting of a 163 
pointing WSRT coverage of M31 (|Braun et al.j2009|) and 



Very Large Array (VLA) coverage of 1V133 consisting of 
6 pointings in the B and CS configurations and a 20 
pointing mosaic in the D configuration. Some early M33 
results base d on the B and CS c onfiguration data were 
presented in Thilker et al. (2002) while an independent 
reduction of all three VLA co nfigurations has been pre- 

The interferometric data 



sented by Gratier et al. 
provides resolution as 



( [2010 ). 

iigh as 15 arcsec in M31 and as 



2007 



of M31 out to 150 kpc ( Ibata et al._ 
SPLASH ch aracterisation of the smooth halo by |Gilbert 
et al.l 2012 1. Furthermore, the data reveals a wealth of 



dwarf galaxies (e.g. Martin et al. 2006 2007 |Richard 



son et al. 


2011), globular clusters and extended clusters 


(Huxor ei 


c al. 


20081 Cockcroft et al. 


2011 Huxor et al. 


2011 Tanvir et al.||2012|); while these extended clusters 



(Huxor et al. 2005) , they appear to have ver y similar 



stellar populations (Mackey et al. 2006 2007|), and do 



no t appe ar to be dominated by dark matter (Collins et 



al. 



2009]. 

he structure around 



high as 5 arcsec in M33. 

2.2.2. H I Properties 

The H I distribution in the extended M31 and M33 en- 
vironment is overlain as contours on the surface density 
of RGB stars in the PAndAS survey in Figure [2j The 
H I contours delineate the diffuse gaseous filament that 
connects M31 and M33 as well as other filamentary fea- 
tures extending both to the northwest of the M31 disk 
and to the southwest. Heliocentric radial velocities of 
various features are indicated with the annotations. The 
intrinsically diffuse nature of these features has been ver- 
ified with follow-up observations using the GBT directed 
at the brightest positio ns along the M31-M33 filament 



M33, appare nt as a distorted 
out er disk, was originally presented by |McConnachie"et 
al. ( 2009 , 2010 1 . This structure is quite clearly asymmet- 
ric about M33, being significantly more prominent in the 
north than in the south. 

The extension of the PAndAS survey towards to north 
of M31, in the region of the two satellite galaxies, 
NGC147 and NGC185, reveals the presence of stellar 
debris that appears to have been torn from NGC147, 



(Braun & Thilker 2004] ) . Despite employing a 25 times 
smaller beam area in these GBT observations, the fea- 
tures were detected at the same low column densities of 
only a few times 10 17 cm -2 . This rules out the possibil- 
ity that the filament is simply a collection of unresolved 
clumps in the discovery observations. Several discrete 
features are also seen in the northern portion of the field. 
It is noteworthy that the M31-M33 filament connects the 
systemic heliocentric velocities of M31 (~ —300 km s _1 ) 
and M33 ( 180 km s" 1 ). 
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but focus ed upon the m a in dis k and inner halo substructure of M31. The left-hand panel is overlain with H I 
e study of Thil ker et al.| j2004[ ) which reveal structure at a spatial resolution of 18 arcmin. The logatithmic 



Fig. 3. — As Figure 
contours drawn from 

contours are drawn at peak brightness temperatures T r = 0.025 to 20 K in steps of 0.15 dex. The right-hand panel is overlain with the 
H I observations of |Braun~e t al. (2009|, further revealing smaller scale structure, at 4 arcmin resolution, in the gaseous component. The 



logarithmic contours are drawn at integrated column densities Njji = 7 X 10 to 4 X 10 cm in steps of 0.28 dex. The annotations on 
both panels indicate the heliocentric radial velocity of H I features, and the image is 108 kpc across. 



The more immediate H I environment of M31 is il- 
lustrated in Figure [3] where the GBT image (left-hand 
panel) and WSRT mosaic (right-hand panel) are simi- 
larly overlain on the PAndAS survey. At these higher an- 
gular resolutions, of 18 and 4 arcmin respectively, much 
of the diffuse gas is not detected but rather only the dis- 
crete High Velocity Cloud (HVC) features within about 
100 kpc of the disk. Of note is that the radial velocity of 
discrete HVC features follows the basic pattern of disk 
rotation with the most negative velocities occurring in 
the southwest and most positive in the north-east. 

The immediate gaseous environment of M33 is illus- 
trated in Figure [4] where the GBT image at 9 arcmin 
resolution (left-hand panel) and the VLA mosaic with 2 
arcmin resolution (right-hand panel) are overlain on the 
the PAndAS survey. M33 is not as rich in discrete HVC 
features as M31 but rather displays a very strong clumpy 
warp of the outer H I disk oriented toward the northwest 
and southeast, strongly suggestive of accretion fuelling 
of the M33 disk by fall-back from the M31-M33 filament. 

3. CORRELATIONS 
3.1. Spatial Correlations 
3.1.1. The disk & halo of M31 

Firstly, we consider the large scale distribution of RGB 
stars and H I gas presented in Figure [2j focusing upon 
material in the vicinity of M31. Leaving aside the inner 
disk region, which will be examined in detail shortly, it is 
useful to consider material along the NW axis, connect- 
ing M31 and M33, and that which lies off this axis in the 
the halo of M31. 

It is clear that that the eastern portion of the halo pos- 



sesses one substantial stellar substructure, namely the 
Eas tern Cloud, as well as a number of dwarf galaxies 
(see Richardson et al.]|2011 ) . This section of the halo is 
effectively devoid of a substantial quantity of H I mate- 
rial, except for a lone clump located to the far north, and 
some lower level material closer to the disk of M31, with 
no significant enhancement in the vicinity of any of the 
dwarf galaxies. Hence, there appears to be no substan- 
tive correlation between the stellar and gaseous material 
in this area of the halo. 

The situation is very similar in the western halo, with 
a one extensive stellar substructure, the SW Cloud, and 
a population of dwarf galaxies, and some potential, but 
low lev el, substructure features. Note that |Mackey et ah] 
(20101 also identified a potential over-density of globular 
clusters, the NW Group, which is not associated with an 
over-density in the stellar density (see Figure [I]). Again, 
this area is almost devoid of significant H I detections, 
other than a large spur at ~ 10° to the major axis, ex- 
tending ~ 65 kpc from the center of M31. Again, there 
is apparently little stellar substructure along the length 
of this gaseous spur, except at its tip where it overlies a 
portion of the SW Cloud. The physical implications of 
this association will be discussed in Section |U 

The most conspicuous stellar and gaseous substruc- 
tures lie along the NW-axis, along the line joining M31 
and M33. In the region south of M31, there is signif- 
icant stellar substructure in the form of the GSS, and 
perpendicular steams, B, C & D, close to the outer disk 
of M31, and the distant stream A. The gas in this re- 
gion is significantly more extended, lying all along the 
axis. There is no significant overlap of H I with the 
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Fig. 4. — The RGB stellar map of M33 and its associated tidal structure, overlain with the H I observations of |Thilker et al.] l |2004| l with 
the GBT, with a resolution of 9 arcmin (left-hand panel), with logarithmic contours are drawn at peak brightness t emperatures Th = 0.15 
to 20 K in steps of 0.18 dex. The right-hand panel presents the VLA observations of H I by l |Thilker et al.|2002| at 2 arcmin resolution. 
The logarithmic contours are drawn at peak brightness temperatures Tg = 0.9 to 30 K in steps ot V.'Zb dex. Tne annotations indicate the 
heliocentric radial velocity of H I features, and the left-hand image is 60 kpc across, whereas the right-hand image is 35 kpc across. 



GSS, with the gaseous material extending from the disk 
of M31 over the inner stream fields and towards M33. 
Significant H I also extends northwards from the disk of 
M31, over the region encompassing the NW stream, al- 
though the gaseous material is more extensively spread 
than the stellar substructure in this area. 

The two panels in Figure [3] zoom in to the outer disk 
region of M31, overlain withH I at the two resolutions 
discussed earlier. While the GSS and Stream C are ap- 
parent as stellar substructures, with Stream B stretch- 
ing back towards M31, there appears to be no significant 
correlation between the stellar material, and H I gas. In- 
triguingly, there is some H I emission, apparent as pop- 
ulation of distinct clouds, aligned with the Giant Stellar 
Stream, but o ffset by ~ 15 kpc; as will be discussed in 
Section |3Xi] this emis sion is also kinema tically corre- 
lated witFthe GSS (see |Braun et al |)2009| ). 

Clearly, the H I associated with the stellar disk pos- 
sesses similar distortions to the underlying stars, includ- 
ing a Northern Spur, located in the outer disk of M31 
(close t o the NE Structure in F igure [I] see also Fig- 



ure 1 in Richardson et al. (2008)) and a potential over 



density hear the Gl-clump; however, it is difficult to dis- 
cern whether either the stellar or gaseous material has 
recently been accreted or is actually disk mate rial that 



ha s been tidally distorted due to interactions (Ibata et 
al] [2005| [Faria et al.|2007 |. 



3.1.2. The disk & halo of M33 

At t he resolution provided by the |Braun fc Thilker| 
( |2004| Figure [2]) , the H I emission from M33 is seen to 
be quite extensive and orientated E-W across the galaxy. 
The stellar substructure, on the other hand, lies mainly 
on the NW-axis, joining the disk of M33 at the northern 
tip of the galaxy. 



The situation becomes clearer in Figure [4] which fo- 
cuses upon the stellar distribution of M33 overlain with 
H I emission seen with t he GBT (left-hand pane l) and 
VLA (right-hand panel) ( |Thilker et al.||2004| |2002| . The 
striking feature in both the H 1 and stellar material is the 
substantial stream pointing towards the NW, as well as a 
less pronounced compone nt in the SE; this ext ended H I 
emission was discussed by Putman et al. (2009), with the 
conclusion that it results from the tidal interaction be- 
tween M33 and M31. A similar conclusio n with regards 
to the stella r component was reached by McConnachie 



et al.| (2009), also presenting a numerical model for the 
tidal interaction of the two galaxies which results in a 
consistent stellar feature. However, while the stellar and 
gaseous material possess a similar alignment, they are 
distinct spatially, with the majority of the gas streaming 
from the Western edge of M33, whereas the stellar ma- 
terial streams from the Northern tip of the galaxy. Fur- 
thermore, while both the Northern and Southern stellar 
and gaseous components are asymmetrical, it is appar- 
ent that the stellar component possesses a much more 
pronounced asymmetry. 

3.2. Kinematic Correlations 

While the radio observations provide a global picture of 
velocities, the nature of multi-object optical spectroscopy 
ensures that the stellar kinematics arc determined in a 
series of discrete fields. As well as the parallel kinematic 
survey undertake n by members of the PAndAS team (e.g. 
Ibata et al.||2004[ ), significant effort has been undertaken 
by other "groups (e.g. the SPLASH survey Kalirai et al 



2009 Gilbert et al. 



2009 Kalirai et al. 20 I'D). However 



the gene ral lack of distinct spatial correlations (see Sec 



tion 



3.1) means that direct comparison of kinematics at 



specific locations is not possible, and only a general com- 
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parison can be undertaken. 

3.2.1. The disk & halo of M31 

The most prominent stellar feature in the halo of 
M31 is the GSS, lying close to the (3-D) distance 
where it meets the spir al disk, sweeping backwards to 
over 



100 kpc behind (Ibata et al.||2001 



McConnachie 



et al.| |2003[ |Fardal et al.||2012p . A~steilar kinematic 



survey oTthe GS S shows a strong velocity gradient 



Ibata et al. 2004); the most distant tip of the stream 
.5 degrees away from M31's disk, have velocities of 
v i^\^ ~ 300 km s _1 , placing it essentially at rest with 
regards to M31, while those approaching the disk are 
travelling at —500 km s _1 («M3i~ —200 km s _1 ). 
While there is no ap pare nt gas associated with the extent 
of the GSS (Section 3.1 I, there is significant parallel H I 
emission. This gas possesses a strong velocity gradient, 
but different to that of the Giant Stellar Stream (Fig- 
ure^; close to the M31 disk, the gas is moving close to 
systemic velocity of M31, whereas away from the disk, 
the velocity becomes more positive. The velocity gra- 
dient bridges the systemic velocities of M31 and M33, 
demonstrating a direct connection and common origin. 

North of M31, along the M31-M33 axis, there is sig- 
nificant H I emission which seemingly overlaps with 
North- Western Stream. The velocity of the H I 
(vh~ —380 km s _1 ) appears to be a continuation of the 
gaseous stream connecting the two galaxies. Presently, 
there are no velocities of the stellar structure in this re- 
gion, and so it is difficult to comment on any putative 
connection, although the nu merical mode l s suggest that 
these featu res are unrelated (Bekki 2008 McConnachie 
et al.|2009 b. 

Figure [3] reveals the kinematic structure of the H I 
in the outer disk and inner halo of M31, where is an 
enhancement of H I close to, but not completely aligne d 
where the GSS meets the disk of M31 (Section 3.1.1). 
Intriguingly, while offset spatially, this gas is moving at 
Vh~ —510 km s _1 , similar to the velocity of the stellar 
content of the GSS, and with «M3i~ —200 km s _1 . At 
higher resolution, this gas is decomposed into distinct 
clumps; this will be discussed in Section [4j 

Also seen in stellar maps underlying Figures [3] are 
Streams C and D. While there is apparently no emission 
associated with Stream C, there is one blob of emission 
almost overlaying Stream D with Vh ~ —400 km s _1 . 
This emission arises close to two kinematic fields on 



Stream D obtained by Chapman et al. (20081, both of 



which yielded velocities of —391 km s _i ; given the sim- 
ilarity of these velocities, and their distinct difference 
from the large scale H I gradient between M31 and M33 
(with an expected velocity of ~ —250 km s _1 ) this is 
strongly suggestive of an association between the stellar 
and gaseous material in Stream D. 

3.2.2. The disk & halo of M33 

There is a gross alignment of the stellar debris and that 
of the H I emission in the outskirts of M33, although 
there again is the s patial offset between the two compo- 
Figure |4j) . The kinematic proper- 



nents (Section 3.1 



ties of the gaseous material within the disk of M33 can 

s In the following, refers to heliocentric velocities, whereas 
v M3l is relative to M31 



be de scribed with a galactic rotation model (i.e. Corbelli 
2003 1 , and also to the emission beyond edge of the stellar 
disk. Unfortunately, the only stellar kinematics obtained 
beyond the opti cal disk were obtained alon g the South- 
ern major axis (McConnachi e et al.||2006), away from 



the p rominent stellar debris ( McConnachie et al. 2009 
20101, and hence no detailed comparison can be made. 



4. INTERPRETATION 

What does the general lack of spatial and kinematic 
correlations between stellar substructure and H I tell 
us about the ongoing accretion in the M31-M33 sys- 
tem? Clearly, the axis connecting these two galaxies 
contains the majority of the substructure, but it ap- 
pears that, in this case, the stellar and H I are the 
result of distinctly different accretion events; the main 
gaseous bridge through this region arises from the in- 
teraction between M31 and M33 alone, with any stellar 
material too sparsely distribute d to be detected (Bekki 
2008| |McConnachie et al.|]2009| ). Given M33's relatively 
large mass, and potentially large impact-parameter, it 
has been able to retain the bulk of its stars and gas dur- 
ing the interaction with M31, although the interaction 
has distorted both the stellar and gaseous disk of M33. 

The progenitor of the Giant Stellar Stream, however 
would have been a lower mass system that has been in- 
teracting with M31, and, due to this interaction any pro- 
genitor gas was lost long ago, while the remnant stellar 
debris exists in the form of tidal streams and shell-like 
caustics. At first glance, the lack of significant gaseous 
structures corresponding to the GSS and other associ- 
ated stellar features may seem surprising. However, as 
noted previously, while the gaseous and stellar compo- 
nents of an orbiting satellite are both subject to gravi- 
tational (tidal) forces, gas is also subject to shocks and 
hydrodynamic drag. One possibility is that ram pres- 
sure stripping led to nearly complete loss gas from the 
GSS progenitor well before the formation of the observed 
GSS and associated structures. However, even if the 
ram pressure is not sufficiently efficient to completely 
strip the progenitor of its gas, hydrodynamic interactions 
can erase potential spatial and kinematic correlations be- 
tween the two components. We are presently studying 
the gaseous and stellar streams using high resolution nu- 
merical simulations (Shannon et al., in preparation) and 
preliminary results suggest that the two components can 
experience significant dislocation over the course of a sin- 
gle orbit, with the loss of angular momentum and energy 
due to hydrodynamic forces causing the gas to sink to 
the center of M31, unlike the much longer- lived stellar 
structures. 

One intriguing correlation occurs where the GSS meets 
the disk of M31 (see Figure [3]). While the gas is spatially 
offset by ~ 15 kpc, it velocity of —500 km s _1 is re- 
markably close to that of the GSS at this point. While 
they may be unrelated, it is worthwhile examining if it is 
plausible that they are physically conected. Such a ques- 
tion is complicated by the fact that there appears to be 
very little gas associated with the GSS, consistent with 



a lack of recent star- formation (Brown et al. 2006), but 



this could be solved if the enhancement ol ti 1 marks the 
(as yet unidentified) progenitor of the GSS, with the off- 
set in the spatial location illustrating the action of ram- 
pressure stripping; such a conclusion would indicate that 
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the accretion of the GSS is a relatively recent event and 
has not undergone several complete orbits, or by a pro- 
genitor massive enough to retain some of its gas through 
a prolonged accretion. This latter option is in agreement 
with the prediction of iFardal et al. ( 2008 ) who suggested 



that the source of the GSS was a rotating disk galaxy 
with a stellar mass of ~ 1O 9 M0. 

An alternate explanation also requires that the progen- 
itor of the GSS was gas-rich, but that some gas remains 
orbiting with the stars. This gas would be too tenu- 
ous to be visible in the observations presented in this 
paper, but as the stream is funnelled into the central 
regions, orbits converge and the local densities of both 
stars and gas increases. Given this density enhancement, 
the gas could become visible to our observations. Again, 
the offset in spatial location probably indicates ongoing 
ram-pressure, although there may be a delicate balance 
between the amount of ongoing ram-pressure and the dy- 
namical shepherding by the associated stream of stripped 
dark matter, so that any can remain associated with the 
stream. Detailed simulations of the formation and evo- 
lution of the GSS, considering realistic gas physics, are 
required to address these questions. 

Away from M31, there is a tantalising correlation be- 
tween H I and stellar material in the vicinity of the SW 
Cloud. There are two potential interpretations for this, 
one with the SW Cloud representing an outwardly mov- 
ing agglomeration of shredded stars, left over from a close 
interaction with M31 ~ 250 million years ago. This in- 
teraction would have to have been with a relatively pris- 
tine, gas-rich progenitor to leave the shredded debris with 
gas. On the other hand, the SW Cloud could represent 
the enhancement at the turning point of a stellar stream, 
and this stream would again have to be gas-rich to ex- 
hibit a similar enhancement in the H I. Both scenarios 
require a relatively recent first passage of M31, ensuring 
that the progenitor does not lose all its gas through ram- 
pressure stripping through an extended interaction. The 
estimated stellar and gas mass in the region of the SW 
Cloud are ~ 10 7 M Q and ~ 5 x 1O 5 M respectively, sug- 
gesting that, if this is representative of a recent accretion 
event, then the progenitor w ould have properties simila r 
to the dwarf irrcg ular IC10 HMcConnachie et al.||2012a[ ). 

While we do not have kinematics for the stellar com- 
ponent of the SW Cloud, it does appear to be spa- 



tially correl ated with three globular clusters (see Mackey 
et al.|[2010| for which velocities have been recently de 



termmed (these will b e presented in detail in f orth- 



coming public a tions by Mackey et al. (2012) and Vel 
janoski et al. (2012)). If we take their velocities, ot 

433 km s" 1 (PAndAS-7), 411 km s" 1 (PAndAS- 

8) and 363 km s" 1 (PAndAS-14), each with an er- 
ror of less than 10 km s -1 , as being representative to 
that of the SW Cloud in general, we can compare this 
gaseous spur velocity of ~ —470 km s _1 ; which this is 
broadly consistent with the globular cluster velocity, es- 
pecially if we expect a velocity gradient along any pu- 
tative stream connecting the two component. It should 
be noted, however, the width of the H I in the SW Spur 
is ~ 200 km s _1 , again indicative that we may be look- 
ing both outbound and inbound gases streams, although 
until we have stellar kinematics for the SW Cloud, and 
detailed dynamical modelling, the link between the SW 
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Fig. 5. — The heliocentric velocity field from a distribution of 
likely M33 orbits overlaid on an image from a simulation of a hy- 
pothetical interaction between M33 and M31. Also shown are the 
positions of some of the main gas features from Fig. [2] and their 
observed line of sight velocities. The simulation contains only stars 
and dark matter without gas and shows an enhanced view of the 
distorted outer disk of M33 and tidally stripped stars. Some outer 
disk stars are thrown outward on a tidal tail to the SE while oth- 
ers form a tidal bridge back towards M31. The velocity field is 
generated from a distribution of orbits generated from a Baycsian 
analysis that are statistically consistent with the known distances, 
radial velocities and proper motions of M31 and M33 and their 
given errors. The orbits are further constrained by the assumption 
that the M33 orbit passed within 30-60 kpc at pericentre assuming 
a mass model of M31 with M = 2.5 X 1O 12 M and r 2 oo = 280 
kpc with the plane of the orbit passing close to the "gas bridge" 
between M31 and M33. The close agreement between the velocity 
of density peaks on the gas bridge and orbital velocities of M33 in 
interacting models support the idea that the bridge may be tidally 
stripped gas from M33 from a recent interaction. A caveat is that 
the pre-encounter gas disk of M33 would need to extend to ~ 15 
kpc. 



Cloud and SW Spur remains circumstantial. 

The distribution of gas within the M31-M33 system, 
especially the warped gas disk of M33 and the appar- 
ent gas bridge connecting M31 to M33, strongly suggests 
a past interaction of the two galaxies. It has been ar- 
gued that the misalignment between the inner and outer 
gaseous disk of M33 and its interpretation as a warped 
disk could be the result of a recent tidal interaction (this 
is discussed in more detail below). The faint extension 
of the stellar disk with a similar alignment as the outer 



gas disk can be repro duced in interaction models (Mc 
Connachie et al.||2009| ). The gas bridge between the two 



galaxies has als o been mode lled in the context of a past 
interaction (e.g.jBckki 2008J . 

While the interaction hypothesis is a plausible explana- 
tion for the phenomenology of the gas in the M31-M33 
system it ultimately depends on knowing the orbit of 
M33. R ecent measurements of b oth the pr oper motion 
of M 33 ( |Brunthaler et al.||2005[ ) and M31 ( |Sohn et al. 
2012 1 provide strong constraints on the M33 orbit within 



assumed mass models of M31. We have be en recently ex- 
tending the dynamical models described in |Mc Connachie 
et al. ( 2009 ) to describe the distorted outer stellar disk of 
1V133 to a range of mass models for M31 and M33 using a 
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Bayesian analysis to determine the distribution of orbits 
that are statistically consistent with the observed dis- 
tance, radial velocities and proper motions of the galax- 
ies. This is a work in progress but we present some pre- 
liminary results here relevant to the observed gas dis- 
tribution. In Figure [5j we present the distribution of 
orbits and the resulting line-of-sight velocity field consis- 
tent with the observations with the constraints that the 
pericentre for the M33 orbit is in the range of 30 — 60 kpc 
in a mass model of M31 with M = 2.5 x 10 12 M Q with 
r 200 = 280 kpc. The orbits are computed within the 
model potentials including a Chandrasekhar drag term 
to model dynamical friction. The drag coefficients are 
calibrated against live N-body simulations to ensure their 
applicability to the Bayesian analysis. We further con- 
strain the distribution to lie within 20 degrees of the 
plane containing the gas bridge with the implicit assump- 
tion that the bridge is created from tidally stripped gas 
from M33. We have also overlaid the stellar distribu- 
tion from one N-body simulation including dark matter 
and stars without gas to orient the orbital distribution 
as well as illustrate the expected distortion of the outer 
disk of M33. We note that the complex stellar features 
seen around M33 in the simulation would be undetectable 
using current observations but we present an enhanced 
view to illustrate the complex dynamics. This model suc- 
ceeds in two points in reference to the gas distribution. 
First, the model predicts that the outer disk of M33 will 
be warped and twisted as the result of an interaction 
while the inner bright disk will remain unscathed. Sec- 
ond, tidally stripped stars result in a tidal tail extending 
to the SE and a tidal bridge that falls back onto M31. 
While we have not modelled the gas in this simulation, 
we expect stripped gas to follow approximately ballistic 
trajectories similar to the stars. The heliocentric veloc- 
ity of the distribution of orbits is consistent with the 
observed velocities of -240 and -280 km/s seen in promi- 
nent gas features in the bridge. A tidal bridge of gas 
falling back onto M31 leading the position of M33 on its 
orbit after a past interaction is a plausible description of 
these features. 

Finally, we look more closely at the stellar and gaseous 
material in the outskirts of M33, with distinct features 
visible in both the NW and SE (see Figure [I]). The posi- 
tion angle of the H I feature is offset from that of the disk 
by ~ 30° — 40°, and has been interpreted as a warped 
dis k in projection and modelle d in terms of tilted rings 



loo p-r 



by |Corbelli fc Schneider| ( |1997j) . The stellar features lie 
roughly between the disk and H I position angles. 

The origin of both the H I and stellar features has 
been attributed to tidal interactions with M31 during a 
close encounte r flPutman et al.|20 09 ; McConnachie et al 
|2009[ but see |Keakes fc JNewton||1978| tor an early dis- 
cussion of this possibility), though models of a steady 
state precessing warp in a mo del with a flat tened halo 
have also been considered (e.g. Kuijken| 1991 ). Compet- 
ing with tidal effects is the gravitational field of the M33 
disk and dark halo. In addition, the H I disk experiences 
ram pressure due to its interaction with the gaseous halo 
of M31. We can distinguish between three regions: The 
outermost part of the disk where tidal interactions and 
ram pressure strip material from M33 leaving behind a 
stream, which roughly traces the M33 orbit; an interme- 
diate region, where material is stripped from the M33 




10 I :■ 

R (kpc) 

Fig. 6. — The various components of force (per unit mass) for 
the restoring force of M33 (magenta) and the disk (green), ram 
pressure (red) and the M31 tidal field (blue). At the tidal radius, 
Rt ~ 20 kpc, the forces due to the M31 and M33 tidal fields are 
equal, while ram-pressure stripping can disrupt the gaseous disk at 
R > 7 kpc and completely strip material from beyond R > 14 kpc. 



disk but remains bound to the galaxy; and a close-in re- 
gion, where the disk is warped but material remains on 
roughly circular orbits. 

In this simplified argument, we consider an element of 
the M33 gas disk at galactocentric radius R with surface 
area Afl and surface density E gas - The tidal force on 
this element due to M31 is 



tid 



(1) 



where R p is the perigalactic distance between M31 and 
M33. In deriving this expression we have assumed that 
R p ^> R and that M31 has a flat rotation curve with a 
circular velocity Vm3i- The ram pressure force is 



(2) 



where p gas is the density of the M31 gaseous halo. The 
gravitational restoring force due to the M33 disk is 



47rG£ disk XL as Af2 



(3) 

where Edisk is the total surface density of the disk. 
If we assume an exponential disk, then Sdisk = 
{Md/^GRd) exp (— R/Rd). The restoring force due to 
M33 as a whole is 



-Ft, re 



V 2 

^E Kas AQ 
R 



(4) 



where we assume that M33 also has a flat rotation curve. 
Setting i^tid = -Ft, res yields an estimate for the tidal ra- 
dius: R t ~ R p (V M 33/Vm3i)- 

The different forces are shown in Figure [6j and for 
purely illustrative purposes, we set Vm3i — 250km s~ , 



V M 33 = 100km s- 1 , V Tcl 
M A = 2.6 x 10 9 M o , R d 



400 km s, R n 



50 kpc, 
= 8 x 



1.1 kpc, and p g 
le to the M31 
and due to M33 itself are equal at the tidal radius 



10 5 cm 3 m,H- The forces due to the M31 tidal field 
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R t ~ 20kpc. Ram pressure can disrupt the disk at 
R > 7kpc and strip material from M33 R > 14kpc. 

While these arguments represent order-of-magnitude 
estimates, they appear to roughly correspond to the ob- 
served structures in M33. However, a full exploration of 
the M31-M33 interaction, especially with regards to the 
stripping of M33, using models including gas, stars and 
dark matter, are required to fully understand the com- 
plex dynamics on display; this will form the basis of an 
upcoming paper (Dubinski et al. 2012). 

5. CONCLUSIONS 

We have presented the spatial and kinematic correla- 
tions of the stellar and gaseous substructure within the 
haloes of M31 and M33. While significant substructure 
is apparent in each component, and there is a gross align- 
ment of significant features, mainly along the axis con- 
necting M31 and M33, there is a distinct lack of correla- 
tion between the detailed structures. 

The lack of a truly global kinematic survey of the stel- 
lar substructure within the halos of M31 and M33 limits 
a detailed comparison of stars and gaseous material, al- 
though the growing number of "key-hole" observations 
with 10-m class telescopes is addressing this. However, 
this situation will be resolved with the advent of wide- 
field spectroscopy on large telescopes, such as the pro- 
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posal to build the ngCFHTQ ushering in a new era in 
understanding galactic archaeology. 

The resulting conclusion of this study, therefore, is that 
the gaseous and stellar substructure within the halos of 
M31 and M33 have been built through a number of dis- 
tinct accretion events, but differing physical processes 
have driven the disruption of the stellar and gaseous 
components of any particular substructure, with the lat- 
ter suffering the additional forces of shocking and ram- 
pressure stripping. Given that, through the use of high- 
resolution numerical simulations, our understanding of 
the differing processes driving the disruption of stellar 
and gaseous material, a comparison of their distribution 
through the halo should allow us to dynamically date the 
accretion; this will be the subject of a future contribu- 
tion. 
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